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ABSTRACT 

The hot accretion flow (such as advection-dominated accretion flow) is usually optically thin in the radial 
direction, therefore the photons produced at one radius can travel for a long distance without being absorbed. 
These photons thus can heat or cool electrons at other radii via Compton scattering. This effect has been 
ignored in most previous works on hot accretion flows and is the focus of this paper. If the mass accretion 
rate is described by M = Mq(t / V out ) - 3 and r out = 10 4 r s , we find that the Compton scattering will play a cooling 
and heating role at r < 5 x 10 3 r s and r >5 x 10 3 r s , respectively. Specifically, when Mq > 0.1LEdd/c 2 , the 
Compton cooling rate is larger than the local viscous heating rate at certain radius; therefore the cooling effect 
is important. When Mq > 2L El ^/c 2 , the heating effect at r out is important. We can obtain the self-consistent 
steady solution with the global Compton effect included only if Mo < ^Edd/c 2 for r out = 50r s , which corresponds 
to L < 0.02LEdd- Above this rate the Compton cooling is so strong at the inner region that hot solutions can 
not exist. On the other hand, for r out = 10 5 r s , we can only get the self-consistent solution when Mo < ^Edd/c 2 
and L < O.OlLEdd- The value of this critical accretion rate is anti-correlated with the value of r out . Above this 
accretion rate the equilibrium temperature of electrons at r out is higher than the virial temperature as a result of 
strong Compton heating, so the accretion is suppressed. In this case the activity of the black hole will likely 
"oscillate" between an active and an inactive phases, with the oscillation timescale being the radiative timescale 
of the gas at r out . 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — quasars: general — 
X-rays: general 



1. INTRODUCTION 

Compton scattering between photons and electrons is an 
important process in astrophysics. If the photons are not pro- 
duced at the same place where the electrons are located, we 
call it "global Compton scattering". Momentum and energy 
of photons and electrons can be exchanged in this process and 
these two aspects often play an important role in determin- 
ing the dynamics of the gas flow. On the galactic scale, this 
so-called radiative feedback mechanism now is believed to be 
crucial for understanding AGN feedback on galaxy formation 
and evolution (e.g., Ciotti & Ostriker 2001, 2007; Murray et 
al. 2005; Hopkins et al. 2005). On a smaller scale, the effect 
of Compton scattering on the dynamics of gas flows surround- 
ing a strongly radiating quasar has been investigated and out- 
flow is produced as a consequence of Compton heating (e.g., 
Proga, Ostriker & Kurosawa 2008). Following the earlier 
work of Krolik et al. (1981), Mathews & Ferland (1987) con- 
sidered the Compton heating effect for the broad-line region 
of quasar. This effect is also important for the standard thin 
disk, if the disk is warped or irradiated by a source above the 
disk plane (e.g., Shakura & Sunyaev 1973; Begelman, McKee 
& Shields 1983; Dubus et al. 1999). If the accretion flow is 
geometrically thick and optically thin, the photons can travel 
a large distance without being absorbed, therefore the global 
Compton scattering effect is in principle important. This is the 
case for spherical accretion and hot accretion flows. The lat- 
ter includes the advection-dominated accretion flow (ADAF; 
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Narayan & Yi 1994; 1995) and luminous hot accretion flow 
(LHAF; Yuan 2001; 2003), two types of hot accretion flows 
corresponding to low and high accretion rates, respectively. 

For spherical accretion, the interaction of momentum be- 
tween photons and electrons sets up a largest possible lumi- 
nosity the accretion flow could reach, namely the Eddington 
luminosity LEdd (but this limit does not apply when the accre- 
tion flow has a non-zero angular momentum; see, e.g., Ohsuga 
& Mineshige 2007). The effect of the energy interaction be- 
tween photons and electrons in a spherical accretion flow has 
been investigated by Ostriker et al. (1976). It was found that 
when the luminosity is larger than a certain value the outward 
energetic photons could heat gas flow significantly so that the 
local sound speed is larger than the escape speed, or, in other 
words, the temperature is higher than the virial temperature, 
thus the accretion is suppressed; and this effect due to energy 
input from the outgoing radiation field occurs for much lower 
luminosities than the momentum (Eddington) limit. 

In almost all of the previous work on the dynamics of hot 
accretion flows, only the "local" Compton scattering effect 
has been considered while the global Compton effect has been 
neglected. Here "local" means that photons are produced at 
the same region with where the electrons locate. This local 
Compton scattering serves as the main cooling mechanism 
of electrons (Compton cooling) and the main mechanism of 
producing X-ray emission (thermal Comptonization). The 
only works to our knowledge considering the global Compton 
scattering are Esin (1997) and Park & Ostriker (1999; 2001; 
2007). Esin (1997) deals with a one-dimensional ADAF and 
finds that the global Compton heating/cooling is not important 
and can be neglected. Park & Ostriker (2001; 2007) deal with 
a two-dimensional flow and focus on the possible production 
of outflow in the polar region because of strong Compton 
heating there. Their conclusion is that Compton heating ef- 
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feet is important in many cases, which is quite different from 
Esin's result. They do not attempt to obtain the self-consistent 
solutions. 

All their work are based on the self-similar solution of 
ADAF (Narayan & Yi 1994; 1995). The discrepancy between 
Esin and Park & Ostriker is likely due to the additional but 
different assumptions adopted. While the self-similar approx- 
imation is quite successful in catching the main spirit of an 
ADAF, it is not a good approximation when we want to calcu- 
late the radiation since order of magnitude error could be pro- 
duced. This is because the self-similar approximation breaks 
down at the inner region of the ADAF where most of the ra- 
diation comes from. When we consider the effect of global 
Compton heating/cooling, obviously it is crucial to calculate 
the exact spectrum from the exact global solution of the accre- 
tion flow. This is the main motivation of the present paper. In 
addition, important theoretical progress on ADAF solutions 
have been made since its discovery. The two most impor- 
tant ones are the presence of outflow (e.g., Stone & Pringle 
2001) and significant direct electron heating by turbulent dis- 
sipation (e.g., Quataert & Gruzinov 1999), which is much 
stronger than the heating by Coulomb collisions between ions 
and electrons. Both of them are important in determining the 
dynamics of ADAFs while they have not been properly taken 
into account in the above works. In the present paper we will 
focus on the one-dimensional case, given that the uncertainty 
of our understanding to two-dimensional structure of accre- 
tion flow is still large. Of course, a full understanding to the 
two-dimensional case is obviously important and should be a 
topic of future research. 

2. THE IMPORTANCE OF COMPTON HEATING OR COOLING IN 
HOT ACCRETION FLOWS 

2.1. Calculation Method 

Consider a canonical hot accretion flow without taking 
into account the global Compton heating or cooling. Out- 
flow is taken into account by adopting the following radius- 
dependent mass accretion rate (e.g., Blandford & Begelman 
1999): 

M = -4-irrHpv = M ( —) , (1) 

V r out J 

where Mq is the mass accretion rate at the outer boundary r out . 
The value of index s describes the strength of the outflow and 
we use s = 0.3 from the detailed modeling of Sgr A*, the su- 
permassive black hole in our Galactic center (Yuan, Quataert 
& Narayan 2003). The energy equations for ions and elec- 
trons are 

/ de e p e dp\ . _ 
PV \~d7~~p I ~d? ) = q +qie ~ q ' (3) 

where £ e j is the internal energy of electrons and ions per unit 
mass of the gas, q K is the Coulomb energy exchange rate be- 
tween electrons and ions, q~ is the electron cooling rate, in- 
cluding synchrotron and bremsstrahlung emissions and their 
local Comptonization, q + is the net turbulent heating rate, the 
value of 5 describes the fraction of turbulent heating which di- 
rectly heats electrons and we use S = 0.5 again from the mod- 
eling of Sgr A* (see also Sharma et al. 2007). 

We first get the exact global solution of the hot accretion 
flow so that we know all the quantities such as density and 



temperature as a function of radius. This requires us to solve 
the set of equations describing the conservations of mass (eq. 
1), energy (eqs. 2 & 3), and momentum (ref. Yuan, Quataert 
& Narayan 2003). The global solution should satisfy the 
outer boundary condition at the outer boundary r out , the inner 
boundary condition at the horizon, and a sonic point condition 
at the sonic point. To calculate the rate of "global" Compton 
heating/cooling of electrons at a given radius r, we need to 
know the spectrum received at r emitted by the whole flow. 
This requires us to solve the radiative transfer equations along 
the radial direction, which is complicated when scattering is 
important. For simplicity, here we deal with the scattering in 
a simple way and write the received spectrum at r emitted by 
the flow inside of r as, 

Here r is the scattering optical depth from r' to r, r = 
f r t arn e dr', and dL v (r') is the emitted monochromatic lumi- 
nosity from a shell at r' with thickness dr' and height H(r'). It 
includes synchrotron and bremsstrahlung emissions and their 
local Comptonization. We approximate the calculation of the 
unscattered part of dL v (r') by solving the radiative transfer 
along the vertical direction of ADAFs adopting a two-stream 
approximation (Manmoto, Mineshige & Kusunose 1997): 

Ait 2 r- 
dL™{r') = — B v [l -exp(-2^r;)]/Jr' (5) 
V3 

Here B v denotes the Planck spectrum, t* = 
(tt /2) x l 2 K v {Q)H{r') is the optical depth for absorption 
in the vertical direction with k„(0) being the absorption 
coefficient on the equatorial plane. The free-free absorption 
and synchrotron self-absorption are included in this way. 
The strength of the magnetic field in the accretion flow is 
determined by a parameter (3 defined as the ratio of the gas 
pressure to the magnetic pressure and we set (3 = 9. For the 
calculation of the Compton scattered part of dL v {r'), we use 
the approach of Coppi & Blandford (1990; eq. 2.2). The 
integration in eq. (4) begins from the black hole horizon 
r s = 2GM/c 2 . 

The spectrum received at r emitted by the flow outside of r 
is (Park & Ostriker 2007), 

The total spectrum received at r is the sum of F™(r) and 

We assume that the electrons have a Maxwell distribution 
with temperature T e (8 e = kT e /m e c 2 ) and the energy of the 
photon before scattering is e = hv/m e c 2 . Since the electrons 
in the hot accretion flow are relativistic at the innermost region 
and the peak photon energy e > 1, to calculate the average en- 
ergy of a scattered photon, we use the following exact form 
which is valid for any photon energy and electrons tempera- 
ture (Guilbert 1986): 

< ei >= e+ 

^fc -C + sinh " G(e« V*exp (=^) dcj>, 

(7) 

with G(e) = go(e)-gi(e) and the cross-section for scattering 

(8) 
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FIG. 1. — Left panel: The ratio of the Compton heating/cooling rate 
to the turbulent heating rate of electrons in hot accretion flows for M = 
Mo(r/\(ftr s f 3 with three accretion rates Mq = 0.1 MEdd (dashed), 1 MEdd 
(dotted) and 2 MEdd (solid). Right panel: The ratio of the radiation tempera- 
ture # x (see text for definition) and the electron temperature Q e . When 8 X is 
larger (smaller) than 8 e , Compton scattering plays a heating (cooling) role. 

Here K 2 (x) is a modified Bessel function of second order and: 
^)^^( f(f _ 2)+1+f , + ^_|_. (9) 

In the Thompson limit, eqs. (7) & (8) are transformed into the 
familiar form of 

4kT e - e m e c 2 
<ei>=e + e 5-^—, (10) 

m e c l 

and 

a{€,9 e ) = a T . (11) 

The number of scattering in a region of the accretion flow 
with unit width in the radial direction and scattering optical 
depth T es = a(e,6 e )n e is 

N = t„, (12) 

with e and n e are the temperature and number density of elec- 
trons in that region. The Compton heating (cooling) rate in 
that region (with unit radial length) is then 

<?co mp = Jn [Fl n (r) + Fr(r)i e ~ < e ei> dv. (13) 

Note that in the above equation we actually use the moment of 
intensity "J" not "F". Following Park & Ostriker (2007), we 
formally define a "radiation temperature" (or "Compton tem- 
perature"; see also Levich & Sunyaev 1970; Krolik, McKee 
& Tarter 1981): 

x 4m e c 2 J[Fl n (r) + F° M (r)]dv' ( ' 

Under this definition, whether the Compton scattering plays 
a heating or cooling effect roughly depends on whether the 
electron temperature e is larger or smaller than 8 X . In the 
Thompson limit, the Compton heating/cooling rate is exactly 
proportional to (8 X - 9 e ). 



It is important to note here that the radiation temperature 
is obtained from the flux distribution by weighting with the 
factor hv. It physically represents the equilibrium between 
Compton heating and cooling. For a typical quasar spectrum, 
it corresponds to 2 x 10 7 K or several keV where the bulk 
of the radiation is emitted in the UV or, in some cases, IR 
portions of the spectrum (Mathews & Ferland 1987; Sazonov, 
Ostriker & Sunyaev 2004). As we will see below (e.g., ref. 
Fig. 3 (c)), the radiation temperature of an ADAF spectrum is 
much higher because of its different spectrum. 

2.2. Results 

The dominant heating term of electrons in eq. (3) is 
#vi Sj e = 5q + . We compare the rate of global Compton heat- 
ing/cooling with q'vis e and the results are shown in Fig. 1 
for M = M (r/10 4 r s ) 03 with Af = 0.1, 1, and 2M Edd (M Edd = 
Lndd/c 2 ). At large radii, r > 5 x 10 3 r s , Compton scattering 
heats electrons; while at small radii, r<5x 10 3 r s , it cools 
electrons. This is of course because the radiation temperature 
8 X is lower (higher) than the electron temperature e at the 
small (large) radii, as shown by the right panel of Fig. 1. 

We can see from the left panel of Fig. 1 that the Compton 
effect is important when Mq > 0.1M E dd- In this case, its cool- 
ing effect can not be neglected. The corresponding accretion 
rate at the black hole horizon is ~ 10~ 2 M Edd and the corre- 
sponding luminosity is <~ 5 x 10~ 4 L Edd . The "lowest" value 
of Mq above which Compton heating effect is important is a 
function of r out . For r out = 10 4 r s , this value is ~ 2M Edd and 
the corresponding luminosity is <~ 2 x 10~ 2 L Edd . When r out is 
larger, the value of critical Mq is lower. In reality r out usually 
has a largest feasible value. If the ADAF starts out from a 
transition from a standard thin disk, r out equals to the transi- 
tion radius. If on the other hand the accretion flow starts out 
as an ADAF such as in our Galactic center, r out should be de- 
termined by the Bondi radius. Outside of Bondi radius, the 
effect of Compton heating is not so clear, because matching 
an ADAF solution to one with proper boundary condition at 
infinity is an unsolved problem. 

Our result that Compton scattering heats electrons at large 
radii while cools electrons at small radii is qualitatively con- 
sistent with both Esin (1997) and Park & Ostriker (2001; 
2007). However, Esin (1997) found that Compton heating 
rate is always even smaller than the Coulomb collision heating 
rate therefore is negligible. This is different from our results 
and Park & Ostriker (2001; 2007). The reason may comes 
from some over-simplifications and the different (old) ADAF 
model adopted in Esin (1997). Overall, we see that for ex- 
tended solutions (r out > 10 4 r s ), both Compton cooling in the 
inner parts and Compton heating in the outer parts dramati- 
cally alter the solutions when L > 10~ 2 L Edd . 

3. THE SELF-CONSISTENT SOLUTIONS 

The above result indicates that we should take into account 
the effect of global Compton heating/cooling when we cal- 
culate the global solution of the hot accretion flow when M 
is relatively large. This has not been studied in Esin (1997) 
and Park & Ostriker (1999; 2001; 2007). We use an itera- 
tion method to achieve this. We first solve the global solu- 
tion without considering the global Compton effect, calculat- 
ing the rate of Compton heating/cooling at each radius as de- 
scribed above, g comp . We then include this term in the energy 



4 



**** 





10 14 16 16 20 



r/r s log|>(IIz)] 



FIG. 2. — The Compton effect for a stellar mass black hole with black hole mass of M = 10 Mq and mass accretion rate M = (557" J ^Edd- The solid and 

dashed lines are for the solutions before and after the global Compton effect is taken into account, (a) The electron temperature profile, (b) The ratio of the rate 
of Compton heating/cooling and the turbulent heating of electrons, (c) The spectrum of the accretion flows. 



equations of electrons: 

/ de e p e dn\ _ 
PV \~dr~~f?~dr ) + * e "? +< ?com P , (15) 

and calculate the "new" global solution of the accretion flow 
based on this "new" equation. Then we get a new Compton 
heating/cooling rate. If the new rate is not equal to the guessed 
value we replace the guessed value with the new one and re- 
peat this procedure until they are equal. However, we must 
emphasize that the solution obtained by the above approach is 
actually not exactly "self-consistent". We use eq. (13) to cal- 
culate <7 comp . But in eq. (13), only the local Compton scatter- 
ing is considered when calculating F™ + F° at , and the global 
scattering is difficult to include because we don't know the 
spectrum emitted at other radii which again requires to con- 
sider global scattering. This is difficult to deal with even using 
the iteration approach. The best way to solve this problem is 
by Monte Carlo simulation combined with iteration method. 
This is beyond the scape of this paper and will be our next 
work. On the other hand, we believe our result should be a 
good zeroth-order approximation to the real solution. 



Bearing this in mind, Figs. 2&3 show the calculation re- 
sults. Figs. 2 (a-c) are for a stellar mass black hole with black 
hole mass M = 10 M and M = (r/50 r s ) 03 M Edd . Fig. 2(a) 
shows the electron temperature of the global solution without 
(solid line) and with (dashed line) the global Compton scat- 
tering effect included. Because Compton scattering plays a 
cooling role at small radii, we see that the electron tempera- 
ture decreases after the Compton effect is taken into account 
as we expect. 

Fig. 2(b) shows the ratio of Compton heating and local vis- 
cous heating of the electrons, q com p/qvis,e, before (solid line) 
and after (dashed line) the global Compton effect is included. 
It is interesting to note that at r > 4r s , the absolute value of this 
ratio is between 1 -4. Since we typically have q m <C #vis,e> this 
implies that the right-hand-side of eq. (15) is negative in that 
region, i.e., the viscous heating of electrons is smaller than 
its radiative cooling (q~ -q CO m V )- In another words, the en- 
ergy advection of electrons plays a heating role, just like the 
ions in the LHAF solution (Yuan 2001). In the inner region of 
r % 4r s , where most of the radiation comes from, the absolute 
value of <7com P /<?vis,e ~ 0.5. We find in this case (-g comp ) — q~. 
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FIG. 3. — Compton effect for a supermassive black hole, with black hole mass M = 10 s Mq and mass accretion rate M = 0.5 [^-^ 

dashed lines are for the solutions before and after the global Compton effect is taken into account, (a) The electron temperature profile, (b) The ratio of the rate of 
Compton heating/cooling and the turbulent heating of electrons, (c) The spectrum of the accretion flows. The two arrows show the values of the average energy 
and energy-weighted energy of photons emitted by the self-consistent solution. 



So in the innermost region the viscous heating of electrons is 
equal to its radiative cooling (q~— q CO mp)- 

Obviously, after taking into account the global Compton 
cooling, the radiative efficiency will increase for a given ac- 
cretion rate. But this does not mean that the highest luminos- 
ity L max a hot accretion flow can produce will increase. The 
main heating mechanism of electrons are viscous heating and 
compression work (the second term in the right-hand side of 
eq. 15) while the main cooling comes from (q~ — q C omp)- The 
highest accretion rate beyond which a hot solution no longer 
exists is determined by the balance between heating and cool- 
ing. The heating term is roughly proportional to M while 
the cooling term roughly to M 2 since Compton scattering is 
a two-body collision process. This is why a hot accretion so- 
lution has a highest M. Obviously, when the global Compton 
cooling <7 comp is included, the cooling becomes stronger com- 
pared to the case of only local cooling q~, thus the balance 
between cooling and heating will occur at a lower M. Ac- 
tually M Q = L Edd /c 2 as shown in Fig. 2 is almost the highest 
accretion rate at which we can get the self-consistent hot solu- 



tion, which is a factor of 2-3 lower than the highest rate when 
the global Compton effect is not taken into account. When Mq 
is higher, we find that we are not able to get the self-consistent 
solution since the flow will collapse due to the strong radiative 
cooling. The decrease of the highest M results in the decrease 
of Lmax- Our calculation shows that L max decreases by a factor 
of <~ 2, i.e., we now have L max <~ 3%LEdd- 

We now check how different the spectrum produced by the 
self-consistent solution is compared to the spectrum produced 
by the "old" solution. Fig. 2(c) shows the spectra from the 
hot accretion flow before (solid line) and after (dashed line) 
the Compton effect is taken into account. We see that both the 
luminosity and the cutoff energy of the spectrum (i.e., the cor- 
responding frequency of the peak of the spectrum) decrease 
because of the global Compton cooling. This is of course be- 
cause the electron temperature of the self-consistent solution 
decreases compared to the "old" solution. We would like to 
emphasize again that only local seed photons are considered 
when we calculate the spectrum although we do consider the 
global scattering in calculating the dynamics. Our calculation 
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shows that -q com p ~ <f, so we expect that when the global 
Compton scattering is considered the luminosity of the "ex- 
act" self-consistent solution will be <~ 2 times higher than that 
shown by the dashed line in Fig. 2(c). But the slope and the 
cutoff energy will not change because they are irrelevant to 
the amount of seed photons. 

Figs. 3 (a-c) are similar to Figs. 2 (a-c), but are for a su- 
permassive black hole with M = 10 8 M Q and accretion rate 
M = 0.5(r/10V 5 ) 3 M Edd . We see from the figures that the 
electron temperature decreases after the global Compton scat- 
tering effect is taken into account as we expect, because in 
most region the Compton scattering will cool the electrons. 
Correspondingly, the luminosity of the accretion flow also de- 
creases by roughly a factor of 2 and the cutoff energy of the 
spectrum also becomes smaller. 

We have also calculated the average energy of the photons 
emitted by the self-consistent solution shown in Fig. 3 (c), h < 
v > = J L v dvj J {L v /hv)dv, and the corresponding energy 
of the radiation temperature 8 X at r ouU hv x = m e c 2 9 x . The 
results are ~ 1 eV and 100 keV respectively, and they are 
shown by two arrows in Fig. 3 (c). These values are much 
higher than that of a typical quasar spectrum where, e.g., hv x 
is only several keV (Mathews & Ferland 1987; Sazonov et al. 
2004). 

The spectrum shown in Fig. 3 (c) extends to very high en- 
ergy, > MeV. Observationally, the e-folding energy of the av- 
erage power-law X-ray spectrum observed by Ginga, OSSE, 
and EXSOSAT of radio-quiet Seyfert Is is E c = 0.7^ MeV 
(Zdziarski et al. 1995; Gondek et al. 1996), which is consis- 
tent with the model given the (large) error bar. Better data is 
required to constrain the theoretical model. 

As we state in §2.2, Compton heating effect at large radii is 
another obstacle for us to obtain the self-consistent solution. 
For Mo = M E( jd, if r out > 10 5 r s , we find that the Compton heat- 
ing effect around r out is so strong that the equilibrium temper- 
ature of electrons would be higher than the virial value defined 
as 5/2 k T y „ = GMm p /r, which will in turn make the ion tem- 
perature also higher than the virial one due to the Coulomb 
coupling between them. In this case, the gas is unbound thus 
can not be accreted. The corresponding highest luminosity in 
this case is ~ 2%Ledd- This value is similar to that obtained 
by Ostriker et al. (1976) and Park & Ostriker (2001). From 
Figs. 1 & 3(b), we expect that Mq and the critical r out (signed 
as r v i r i a i) beyond which the equilibrium temperature is higher 
than the virial temperature are roughly anti-correlated, i.e., a 
lower Mq corresponds to a larger r^ai- Exact estimation of 
the relation between Mq and r^ai is not straightforward. This 
is because we need to know the radiation temperature 8 X as a 
function of Mq which requires numerical calculations. Note 
that from Figs. 1 & 3(b) the minimum value of r v i r ; a i should 
be larger than ~5x 10 3 r s . 

Although no steady self-consistent solution exists due to 
the strong Compton heating at and beyond r v iriai, an "oscilla- 
tion" of the activity of the black hole is expected (e.g., Cowie, 
Ostriker, & Stark 1978; Ciotti & Ostriker 2007). When the 
accretion rate is high, only the gas inside of r v ; r i a i can be 
accreted. This active phase will last for a timescale of the 
accretion timescale at r v iriai- Then all the gas will be used 
up and the active phase stops. In this case, Compton heat- 
ing also stops so the gas outside of r v ; r i a i will be cooled by 
radiation and be accreted again and the cycle repeats. The 
time the non-active phase will last is determined by the ra- 
diative timescale of the gas at r V mai, since it is longer than 



the accretion timescale there. An alternative consequence of 
the strong Compton heating at large radii is that the accretion 
can be self-regulated by irradiating the outer flow (Shakura & 
Sunyaev 1973). That is, the strong Compton heating will not 
completely stop the accretion, but only decrease the accretion 
rate. This then reduces the energy release in the inner part, 
which in turn reduces the irradiation. A multi-dimensional 
numerical simulation is required to solve this issue and accu- 
rate time-dependence is needed as well since steady solutions 
may not be stable. 

For the massive black holes seen in the nuclei of most 
galaxies the Compton heated interruption of the high luminos- 
ity states should be typical if a hot accretion flow exists there. 
We now know that the accretion flow in low-luminosity AGNs 
is of this type (see Yuan 2007 and Ho 2008 for reviews). For 
luminous AGNs such as quasars, although people incline to 
think it is a standard thin disk which is optically thick, many 
problems remain for this model (e.g., Shlosman, Begelman, & 
Frank 1990; Koratkar & Blaes 1999). If the actual accretion 
flow is radially optically thin to Compton scattering, similar 
to the hot accretion flow, our analysis applies. 

This kind of oscillation does not apply to stellar mass black 
holes in our Galaxy. This is because the prerequisite for such 
oscillation is that the accretion rate is large and the hot accre- 
tion flow extends to large radii. For a stellar mass black hole, 
the accretion material comes from the companion star and it 
starts out as a standard thin disk. In the hard state the stan- 
dard disk does not extend to the innermost stable circular orbit 
but is replaced by a hot accretion flow within a transition ra- 
dius r tr . However, when the accretion rate is high, r u is small 
(Yuan & Narayan 2004). So Compton scattering cools rather 
than heats the hot accretion flow. But on the other hand, the 
photons emitted by the hot accretion flow will heat the cool 
electrons in the standard disk at r tr and thus will change the 
dynamics of the transition. This effect has never been noted 
and could be a topic of future work. 

4. SUMMARY AND DISCUSSION 

For a geometrically thick and optically thin hot accretion 
flow, the photons can travel for a long distance without be- 
ing absorbed, and thus be able to heat or cool electrons via 
Compton scattering. We investigate this global Compton scat- 
tering effect and find that for an accretion rate described by 
M = Mo(r / r out ) 03 the Compton cooling effect will be impor- 
tant when Mq > O.lLedd/c 2 ; while the Compton heating ef- 
fect will be important when Mq > 2L£dd/c 2 and r out = 10 4 r s . 
Specifically, the scattering heats electrons at r > 5 x 10 3 r s 
while cools electrons at r < 5 x 10 3 r s . If r out is larger, the 
critical Mq above which the Compton heating effect is impor- 
tant will become lower. 

We have successfully obtained the self-consistent steady so- 
lution with this effect included for Mo < Z^Edd/c 2 and r out = 
50r s . But when Mo > ^Edd/c 2 and L > 2%Lndd we fail be- 
cause of the strong radiative cooling (local plus global Comp- 
ton scattering). It is also difficult to get the self-consistent 
solution when Mq > L Edd /c 2 (L > l%L Edd ) and r out > 10 5 r s . 
This is because in this case the Compton heating is so strong 
at and beyond r out that the equilibrium electron temperature 
there will be higher than the virial temperature. More gener- 
ally we expect that the radius where the equilibrium temper- 
ature due to the Compton heating is equal to the virial tem- 
perature, ryiriai, is anti-correlated with M out . We argue that the 
black hole will manifest an oscillation of the activity in the 
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case that we fail to get the steady solution. The period will be 
the radiative timescale of the gas at rvinai- 

All our discussions so far are for a one-dimensional (but 
not spherical) accretion flow. Although big uncertainties ex- 
ist for the vertical structure of accretion flow, we are cer- 
tain that when M is high, the scattering will be important, 
and consequently much of the luminosity will "leak out" per- 
pendicular to the accretion flow as in the standard thin disk. 
This will have two effects. One is that the highest luminos- 
ity up to which we can get the self-consistent solution with 
the global Compton effect included will be higher. In addi- 
tion, the Compton heating will be stronger in the vertical di- 
rection than in the equatorial plane of the flow. As a result 
strong wind will be launched as pointed out by Park & Os- 
triker (2001; 2007) and found by Proga et al. (2008). All of 
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the described effects are likely to become significant for AGN 
accretion flows having L > 10~ 2 LEdd and optically thin (ref. §3 
for discussion of this possibility in luminous AGNs), which, 
we know from recent applications of the Soltan argument (Yu 
& Tremaine 2002), are the phases during which most massive 
black hole growth occurs. 
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